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ABSTRACT 

We analyzed a long duration flare observed in a serendipitous XMM-Newton detection of the 
M star CD-39 771 7B (TWA 11B), member of the young stellar association TW Hya (~ 8 Myr). 
Only the rise phase (with a duration of ~ 35 ks) and possibly the flare peak were observed. We 
took advantage of the high count-rate of the X-ray source to carry out a detailed analysis of 
its spectrum during the whole exposure. After a careful analysis, we interpreted the rise phase 
as resulting from the ignition of a first group of loops (event A) which triggered a subsequent 
two-ribbon flare (event B). Event A was analyzed using a single- loop model, while a two-ribbon 
model was applied for event B. Loop semi-lengths of ~ 4i?„ were obtained. Such large structures 
had been previously observed in very young stellar objects (~ 1 — 4 Myr). This is the first time 
that they have been inferred in a slightly more evolved star. The fluorescent iron emission line at 
6.4 keV was detected during event B. Since TWA 11B seems to have no disk, the most plausible 
explanation found for its presence in the X-ray spectrum of this star is collisional- or photo- 
ionization. As far as we are concerned, this is only the third clear detection of Fe photospheric 
fluorescence in stars other than the Sun. 



Subject headings: stars: activity — stars: 
individual (CD-39 7717B, TWA 11B) 

Introduction 



Among all the processes manifested in stellar 
coronae, flares are the most energetic ones. They 
are supposed to be the result of the energy re- 
lease from m agnetic field reconnection in the lower 
corona (e.g. iKopp fc Pneumanlll976l ). As a con- 
sequence, electrons and ions in the reconnection 
region are accelerated downwards, along the mag- 
netic field lines, toward lower atmospheric layers. 



coronae — stars: flare — stars: pre-main sequence — stars: 



When they reach the upper chromosphere, the lo- 
cal gas is heated and evaporated into the new- 
formed magnetic loops. Thus, the density and 
temperature of these loops increase, causing in- 
tense emission of soft (< 10 keV) X-rays. 

Stellar flares have b een observed i n alm ost all 
the H-R diagram (see IVaiana et al.l Il98l[ for a 



review). However, they are more frequent in late- 
type stars, where they present a large variety of 
sizes and durations. In late-K and M dwarfs 
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(the so-called UV Ceti-type stars; see IPettersenl 



1991L for a description of their main properties), 
moderate flares are frequently observed. During 
such events, the X-ra y flux usually increases b y 
a factor of 2 - 4 (e.g. iRobrade fe Schmitd [20051) . 
Giant flares, in which the X-ray flux increases 
from dozens to hundreds times the quiescent state 
value, have been detected in some M dwarfs 



such as EV Lac dFavata et all l2000h. EQ Pee 
(iKatsova et al. 20021 ). and Prox Cen ( Giidel et al 



2004). The duration of those flares are of the 
order of a few kiloseconds, even in the more ener- 
getic ones. For instance, the giant flare observed 
in EV Lac lasted ~ 5 ks. Long-duration flares are 
more common in pre-main sequence stars than 
in the UV Ceti-type ones. For example, dur- 
ing the 13 days observing run of the ORION 
Nebula Complex by Chandra, at least 19 flares 
with durations abo ve half a day were detected 
" 20051) . Long-duration flares were 



(Favata et al 



also observed in the Taurus star-forming com - 
plex (|Franciosini et alj|2007t IStelzer et alj|2007ft . 
The longer duration of these events is usually 
attributed to the presence of larger coronal struc- 
tures in young stars. 

The XMM-Newton and Chandra missions have 
contributed enormously to the understanding 
of the processes involved in the X-ray emis- 
sion of late- type stars. In particular, the im- 
proved temporal and spectral resolution, to- 
gethe r with the development of theoretical models 
(e.g. IKopp fe Polettol 11984 IPoletto et all Il988t 
Serio et al.l I l991t IGiidel et al.l 11999c iReale et al.l 
20041 : lRealell2007[ ). have provided us with pow- 
erful tools for investigating coronal flares. De- 
tailed diagnostics of X-ray fl ares have been car - 
ried out by different authors. iFavata et al.l (|2005l) 
determined general properties of flaring loops, 
in terms of temperature and semi-length, for 
young stellar objects in the Orion Nebula Com- 
plex (ONC). ^From their analysis, the authors 
inferred loop semi-lengths comparable to the stel- 
lar radius in some cases. They speculated that 
these large structures are connected with the 
proto-planetary disk and are, in fact, the same 
structures that channel the plasma producing ac- 
cretion. A similar work was done for the Tau - 
rus star-forming region ( Franciosini et all [2007). 
where the authors remarked that some coronal 
loops extended up to a distance comparable with 



the stellar radius (see also iGiardino et al.l 120041) . 
Studies for more e v olved stars were done b y, e.g . , 
Reale et al.l (l2004h. ICrespo-Chacon et al.l (|2007n . 
and iTesta et all (|2007l ).~ who found flaring loops 
with semi-lengths of ~ 0.2 — 0.5 (i.e., a rela- 
tively compact flaring corona). 

In the works presented in the previous para- 
graph, the diagnostic of X-ray flare s was done 
using the p rocedure described by IReale et"aT 



(|l997l . 120041 ) to model the decay phase, which 
assumes the flare to be produced in a single 
loop where heating does not entirely drive the 
flare decay. Reale (|2007l ) extended this method 



to the rise phase and compared the parameters 
determined in this way with those obtained by 
analyzing the decay phase for three stellar flar- 



ing loops (see also IPan et aT 1997 ). finding a 



good agreement between them. On the other 
hand, the solar t wo-rib bon model developed by 
IKopp fc Polettol (Il984l) . or its stellar version 
(jPoletto et al.l 1988 ). should be used when heat- 
ing t otally drives the flare evolution ( Reald 2002. 
20031 ). In this model, the reconnection energy 



is supposed to be dissipated immedia tely after 
being released. Kopp fc Polettol ( 19841) assumed 
that only a fraction of the magnetic energy re- 
leased by the reconnection process is used to sup- 
ply th e thermal ene r gy of the newly formed flare 
loops. IPoletto et al. (1988) assumed that a factor 
of 10 % of this thermal energy escapes into the 
X-ray regim e, as suggested by de tailed studies of 
solar flares (ICanfield et al.l ll980T) We also refer 
the reader to IGiidel et al.l (| 1999T ) . who included 
time-dependent conductive and radiative losses in 
the X-ray corona self-consistently. 

In this work, we analyze an XMM-Newton 
serendipitous observation of the young M-type 
star TWA 11B in which the X-ray emission suf- 
fered a continuous increase during approximately 
35 ks. The duration and statistics of the observed 
rise have allowed us: (i) to carry out a detailed 
spectral time-analysis; and (ii) to derive proper- 
ties of the star's magnetic configuration by using 
both the single-loop and the two-ribbon flare mod- 
els described above. 

2. Observation and data treatment 

The XMM-Newton observation (ID 0006220201) 
was perfomed in the revolution 197, between 
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Fig. 1. — Chandra HETG serendipitous detec- 
tion of TWA 11B. The position of the optical 
counterparts of TWA 11 A and TWA 11B arc 
marked. Contours of the PSF observed in the 
XMM-Newton EPIC exposure are overplotted. 



2001 January 4 and 5, for a total duration of 
43 ks. The field is centred on the coordinates 
a 2 ooo = 12 h 35 m 34 s and <5 2 ooo = -39°54'55" (the 
main target being the Seyfert 2 galaxy NGC 4507). 
The EPIC cameras were operated in Full Frame 
mode using the Thick filter. Our target is situ- 
ated 5.18arcmin at the North-East of the central 
source, in the field of view of both the MOS and 
PN detectors. The EPIC-XMM source coincides 
with the position of the optical counterpart of 
TWA 11B. The primary star, TWA 11A (an AO- 
type star), is situated at ~ 8 arcsec from TWA 
11B. Although none X-ray emission is expected 
from an AO star, we investigated the possibility 
that the primary star produced some X-rays that 
could affect our results. We studied a serendip- 
itous detection of HETG Chandra of the source 
(Observation ID 2150). In this 140 ks exposure 
Chandra observation, the X-ray source is clearly 
identified with TWA 11B (see Fig. Q}, while none 
X-ray emission coming from TWA 11A is detected. 
Therefore, we are confident that the X-ray emis- 
sion detected in the XMM-Newton observation 
comes only from the M star TWA 11B. 

The data reduction followed the standard oper- 
ating procedure. We used the version 7.1.0 of the 



XMM-Newton Science Analysis System (SAS) to 
derive a table of calibrated events in the energy 
range 0.3 - 10.0 keV. To extract the events, we 
chose a radius of 30 arcsec, which is slightly larger 
than the 3cr-level of the source's PSF. This assures 
us to lose less than 1% of the counts from the 
source. Different filters were applied to eliminate 
bad events and noise. Note that the observation 
was neither affected by pile-up nor by high flaring 
background periods. 

The X-ray light curve (Fig. _} showed a to- 
tal increase in the star's count-rate of a factor of 
4.2 from the lower level (« 0.25 counts/s) to the 
maximum observed emission. Similar relative in- 
creases in flux were p reviously reported in flare s 
from other stars (e.g. Robrade &: Schmitt 2005). 



A first increase in the light curve is observed 5 ks 
after the beginning of the observation, reaching a 
local maximum only 5 ks later. After a brief de- 
crease, the count-rate continues increasing until 
reaching the global maximum in the light curve 
(this happened 37 ks after the beginning and only 
3 ks before the end of the observation) . The total 
duration of the enhancement is 32 ks (« 9hour_]). 
Note that the observed maximum may or may not 
be the flare peak. In the latter case, the duration 
of the rise phase would be even longer. 

In Fig. _] we plot the evolution of the hardness 
ratio, which is a tracer of the temperature evolu- 
tion, during the observation. Here, the soft energy 
band is defined as the range 0.3-0.8 keV and the 
hard energy band as 0.8 - 4.5 keV. Fig. _] shows 
that the mean coronal temperature reached a max- 
imum at t ~ 7 ks, maintained high with significant 
fluctuations during approximately the next 27 ks 
of exposure, and then it began to decrease gradu- 
ally. 

Overall, the rise phase is unusually long com- 



Pan et al. 1997: Gudel et al.l 1999: 


Reale! 20021 


Robrade & Schmitt 2005; Realc 2007: 


CresDO-Chacon et al. 


2007). To date, so long rise phases have been ob- 



served only in some T Tauri stars (jFavata et al 



1 The rotation pe riod of TWA 11B is not known. 

ISchoIz et al.l (120071 ) measured a projected rotational veloc- 
ity vsini = 12.11 ± 0.93 kms" 1 , what leads to an upper 
limit in the rotational period of 2.7 days. Thus, the ob- 
servation cover, at least, 17% of the rotational period of 
the star. Nevertheless, no indication of occultation of the 
flaring region is observed. 
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2007) to a factor 10 (V892 Tau; iGiardino et al 



Fig. 2.— EPIC-PN light curve of TWA 11B in 
the energy range 0.3 — 10.0 keV. The curve was 
binned to a 1 ks time resolution. The event list 
was corrected of bad events and noise. The expo- 
sure was also corrected of live-time and good-time 
intervals. The gray segments mark the periods in 
which the intrument was turned off. The contin- 
uous line at the bottom is the background light 
curve. The time blocks used in the spectral anal- 
ysis are plotted as dashed vertical lines. 
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Fig. 3. — Hardness ratio evolution of TWA 11B 
during the observation with EPIC-PN. The time 
blocks used in the spectral analysis are plotted as 
in Fig. [2] for clarity. 



20051 : IGiardino et al.l200"llFranciosini et al.l2007l ) 



some of them w ith accretion disks, and several 
RSCVn systems (|Testa et al.l2007tlNordon fc Behar 
2007). In such stars, the relative increase in flux 



2004) with rise-phase times from 8 to 55 ks. TWA 



11B is a weak-line T Tauri M2.5 star with an age 
of ~ 8 Myr (estimated from its membership in 
the TW Hya Association) . While it shows no sig- 
natures of an accretion disk (near-infrared colo r 
excess or strong Ha emission; Stauffer et al.lll995h . 
it is still contracting. Thus, physical conditions in 
its atmosphere should be more similar to those of 
sub-giant stars than of main-sequence ones. This 
may be the reason why long duration rise phase 
flares are observed in both T Tauri stars and RS 
CVn systems. 

A more detailed inspection of the light curve 
shows that the emission increased faster initially, 
it had a well-defined local maximum at f w 10 ks, 
and then it growed again, but more gradually this 
time, until the end of the observation. At the very 
end, there is a hint that the emission was stop- 
ping to increase. Such a long rise phase suggests 
that we observed an uninterrupted sequence of 
flare events involving an extended coronal region. 
However, the faster initial rise and the local emis- 
sion peak - coupled to the earlier hardness ratio 
peak - resembles the evolution of a self-standing 
flaring episode possibly occurring (at least during 
these initial stages) in a single loop (e.g. Realei 



2007), whose decay merged with the ignition of 



the later overlapping events. Fig. |4] demonstrates 
that this peak is indeed significative and cannot 
be attributed to the effect of noise in the light 
curve. This figure shows the cumulative distribu- 
tion of counts for 200 simulated constantly increas- 
ing light curvetH (for t > 5ks) with Poissonian 
noise (shadowed region) together with the cumu- 
lative distribution of counts observed during the 
rise (continuous line). Clearly, the simulations do 
not reproduce the observed distribution of counts 
for the first peak in the light curve. The differ- 
ences are high enough to justify the treatment of 
this event as independent. In the following, we 
work under this hypothesis. 

The large number of counts collected from the 
source in the EPIC cameras allowed us to di- 
vide the observation in several time intervals with 
enough signal to perform a reliable spectral analy- 



ranges from a factor 3 (V410 Tau; lFranciosini et al 



2 Constantly increasing light curves have been chosen for the 
simulations since the global enhancement mostly shows a 
linear pattern. 
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Fig. 4. — Normalized cumulative distribution of 
counts in the light curve of TWA 11B for t > 5ks 
(continuous line). The shadowed region contains 
cumulative curves of 200 simulated constantly in- 
creasing light curves (also for t > 5 ks) with Pois- 
sonian noise. 

sis (> 1000 counts in the quiescent and intervals 1 
and 2, and > 2000 counts in the remaining inter- 
vals after background subtractiorQ) . This permit- 
ted: 1) to investigate the nature of the processes 
taking place in our target; and 2) to derive physi- 
cal properties of it. 

3. Spectral analysis 

For the spectral analysis, the observation was 
split into 10 time-intervals (vertical dashed lines 
in Fig. [2]). These intervals sample different fea- 
tures in the light curve. We used the XSPEC 



spectral fitting package (jArnaudl Il996l |200J) in 
the PN, MOS1, and MOS2 detectors simulta- 
neously. We adopted th e Astrophysical Pla sma 
Emission Code (APEC, ISmith et al.l l2001al) in- 
cluded in the XSPEC software. APEC calcu- 
lates spectral models for hot, optically thin plas- 
mas using the Astrophysical Plasma Emission 



3 Intervals 1 and 2 have been studied separately to investigate 
the evolution of the parameters of the flaring plasma in its 
very first phases. 



Quiescent stote 




Fig. 5.— Observed EPIC PN and MOS spectra 
of TWA 11B during the quiescent state (marked 
with Q in Fig.HJ). 



Database (APED. ISmith et al.lEoOlbh . that con- 
tains the relevant atomic data for calculating both 
the continuum and line emission. Interstellar ab- 
sorption was taken into account using the inter- 
stellar photo-electric absorptio n cross-sections of 
Morrison fc McCammonl (1983), also available in 
XSPEC. 

3.1. Quiescent state 

The lowest (constant) count-rate level, found 
at the beginning of the observation, is assumed 
to be the quiescent state. The X-ray luminosity 
of the star during this lapse of time (assuming 
the same distance than the primary d = 67 ± 3 
pc) is logLx[ergs _1 ] = 29.35, corresponding to 
log(Lx/£boi) = —3.1. These values are typical of 
both pre-main sequence stars members of y oung 
stellar associations (e.g. Kastner et alj|2()()3 i and 



older field M dwarfs showi n g high X-ray activity 
dRobrade fc Schmittl 12005 : Crespo-Chacon et al 
200jj)- The best fit to the X -ray quiescent spec- 



trum (Fig. [3]) is given by a plasma model with two 
equally- weighted thermal components (the output 
parameters of this fit are summarized in Table [TJ . 
Although high uncertainties are obtained for the 
hydrogen column density in our fit, the value of 
iVn obtained by us is very similar to that deter- 
mined for other members of the TW Hya Associa- 
tion (e.g. IStelzer fc Schmittl 12004 lArgiroffi etaL 
2005). Besides, the interstellar extinction was 



found to be negligible in the general direction of 
TW Hya (|Rucinski fc Krautter|[l983l ). 
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Table 1: Output parameters from fitting the X- 
ray EPIC spectra of the quiescent state with a 
2T-model. 



NH 


= 1.7±?-? x 10 20 


cm 


-2 


Z 


= O.17tg:gZ 






kT x 


= 0.27™ keV 






EMi 


= 2.2tH x 10 52 


cm" 


-3 


kT 2 


= 0.98™ keV 






EM 2 


= 2.2±°j x 10 52 


cm" 


-3 


X r 2 cd (d-o.f.) 


= 0.86 (176) 







3.2. Flaring state 

Synthesized stellar-like spectra of solar flares 
(previously subtracted by the quiescent "back- 
ground" spectrum) are generally well-fitted with 
a single therma l component (1-T model, see 



Reale et al.ll2001l ) . Consequently, to analyze phys- 
ical properties of the flaring plasma in our observa- 
tion, we first subtracted the spectrum of the qui- 
escent state from the observed spectrum in each 
time-interval. Then we fitted a lT-model, leav- 
ing its temperature (fcTp) and emission measure 
(EM-p) as free parameters. The values of the col- 
umn density (NH) and abundance (Z) were fixed 
to those derived for the quiescent state. With 
this technique, we interpret that the coronal spec- 
trum during the flare results from adding to the 
quiescent spectrum a third thermal component, 
which is a scribed to heated material filling the flar 



ing lo ops (jReale et al.ll2001t ICrespo-Chacon et al 



2007). Models with additional temperature com- 
ponents were checked, but did not improve the fit 
results significantly. 

Results from spectral fitting are given in Ta- 
ble [2j As expected from the hardness ratio evo- 
lution (Fig. [3]), the temperature peaked in time- 
segment 1. Then, it decreased gradually but al- 
most continuously until the end of the observa- 
tion. On the other hand, the emission measure 
increased continuosly from the beginning of the 
rise and starts to decrease only near the end of 
the observation. 



3.2.1. Evolution of the emission measure distri- 
bution 

In order to approximate physically more re- 
alistic continuous emission measure distribu- 
tions (EMDs) of the plasma^ we used a multi- 



temperature model as in iRobrade fc Schmitt 
(2005). In that work, the authors used a 6T- 
model on a logarithmic, almost equidistant grid 
with temperatures fixed at 0.2, 0.3, 0.6, 1.2, 2.4, 
and 4.8 keV (which correspond to 2.3, 3.5, 7.0, 
14.0, 28.0 and 56.0 MK), sampling those spec- 
tral regions where the XMM-Newton detectors 
are more sensitive. To fit our spectra, the values 
of NH and Z were fixed to those previously de- 
termined for the quiescent state. Thus, the only 
variables are the emission measures. Note that in 
this case we are fitting each spectrum as a whole 
(i.e., the quiescent spectrum was not subtracted 
from the rest of the spectra). Results from apply- 
ing this 6T-model to our observation are shown 
in Table [3] and Fig. [6j and are summarized in the 
following items: 

• The amount of plasma emitting at temper- 
atures above 28 MK is negligible in the qui- 
escent state. 

• Plasma emitting at the highest temperatures 
(~ 56 MK) appears at the beginning of the 
first (faster) rise (time-segment 1), while the 
rest of the EMD curve remains as in the qui- 
escent state. Note that in this time-segment 
we obtained the maximum temperature for 
the flaring component when fitting it with 
the lT-model (see first part of § 13.21 and Ta- 
ble HJ}. 

• The only significant difference between the 
EMD of the quiescent state and that de- 
rived for time-segment 2 (where the light 
curve of the first flare peaks) is an excess 
of plasma emitting at temperatures around 
28 MK. We interpret it as the cooling of the 
material at higher temperatures that was de- 
tected in time-segment 1 together with ad- 
ditional plasma evaporated after the tem- 
perature peak was reached in time-segment 
1 (note that the total emission measure of 
the flaring component in time-segment 2 is 
higher than that measured in time-segment 
1 , as it can also be seen in Table [2]) . 
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Table 2: Spectral results for the flaring component in each time-segment of the light curve (excepting the 
quiescent time-interval). 



Time-segment 


Time- interval 


Central time 


kT F 


EM-p 


X?cd ( d -°- f -) 




(ks) 


(ks) 


(keV) 


(10 52 cm" 3 ) 




1 


5-8 


6.5 


stf 


n 7 9+0.20 
u - ' z -0.17 


1.05 (129) 


2 


8-11 


9.5 


2-4±8:S 


4 o+ - 3 


0.87 (215) 


3 


11 - 17 


14 


o c+0.4 
z -°-0.3 


o OO+0.27 
^•°°-0.28 


1.08 (369) 


4 


17-22 


19.5 


9 97+O.25 
z - z 1 -0.27 


c o+0.3 
°-°-0.3 


0.87 (355) 


5 


22 - 27 


24.5 


1 ss+ - 14 

1.88_ 14 


7 9+O.3 
' ' Z -0.3 


0.89 (406) 


6 


27-32 


29.5 


1 66+ 018 
i - DD -o.o8 


o q+0.3 


1.03 (430) 


7 


32 - 36 


34 


i O-.+0.05 
1 - ol -0.05 


9 2 +0 - 4 

a - z -0.4 


0.97 (379) 


8 


36 - 38 


37 


1 99+O.O6 
i - zz -0.06 




1.01 (243) 


9 


38 - 41 


39.5 


o-99±8:8S 


9-718:1 


0.86 (313) 



Table 3: Output parameters (emission measures, in units of 10 52 cm -3 ) from fitting the whole X-ray EPIC 
spectra to a 6T-model with temperatures fixed at 0.2, 0.3, 0.6, 1.2, 2.4 and 4.8 keV. EM[ is the emission 
measure corresponding to the thermal component of 0.2 keV, EM^ is that of 0.3 keV, and so on. Uncertainties 
are calculated for a 2.7 a confidence level. 

Time Interval EM[ EM!, EM^ EM^ EM^ EM' & 

Quiescent l.OOlgJ? 1.00±l"gg 0.77ljj-{g 1.63±g"|| 0.00t " 00 0.00±° °g 
1 0.66l8: 8 6 8 0.59±S;I| 0.75t ; 7 3 5 l-85l8:g 0.00™ 0.85i° : 29 

9 n1 o+l.24 9 1 ,+0.96 n o-i+1.15 010 +1.15 o 7n +0.66 n nn +1.23 

z u.io_ 13 z.i+_ 214 u.oi_ 31 z.iz_ og7 J-'U_ 2 .i7 u.uu_ 00 

o i 01 +0.61 n 99+1.58 9 CiA+ 0A8 9 7Q+0.54 n ftft +0.61 i Q1 +0.29 

O 1.01_ 1 og U.ZZ_ 22 z - u *-0.78 z -' J -0.60 U ' UU -0.00 1 - yl -0.28 

4 0.85l8:g 1.88l?;» 0.0ll°;» 4.48l°:g l-32±i : g 1.48±i:» 

5 0.54™ 2.42±»:S 0.62™ 5.19™ 1.69li;» 1.41±i:g? 

6 0.92™ 1.94±i;^ 0.96™ 5Al+_\f 2 2.53™ l-01±i:g? 

7 1.26li:« l-36i; : g 1.30li:«| 8.35±1;^ 0.86±J;S 0.58l°;g 

8 2.64™ 0.25™ 3.39™ 8.17™ 0.091^ 0.81±°|? 

9 1.38li:g 0.83™ 4.08li:°S 8-391S.81 0.00™ 0.23t° ;« 



• The excess of plasma emitting in the region 
around 28 MK that was observed in time- 
segment 2 continues cooling towards lower 
temperatures during time-segment 3. At the 
same time, a high quantity of plasma ap- 
pears again at the highest temperatures (~ 
56 MK). This quantity is even larger than 
that measured in time-segment 1 and is ap- 
proximately coincident with the beginning of 
the second (more gradual) rise. 



• During time-segments 4-7, the amount of 
plasma emitting at the highest temperatures 
remains still high, while a large excess of 
material (compared to the quiescent state) 
appears also between 7 and 28 MK (reach- 
ing the maximum always at a temperature 
w 14 MK). The existence of intense, sus- 
tained heating is needed to explain the main- 
tenance of plasma at temperatures around 
56 MK. The evolution of the rest of the EMD 
curve can be interpreted as the cooling of the 
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Fig. 6. — EMDs as derived with the 6T-model. 
The curves are shifted in the y-axis to better ob- 
serve the flare evolution in time. Each one of 
the EMDs (shown as continuous lines) is labeled 
with the same number as its corresponding time- 
segment (see Table [2]). The quiescent is marked 
as Q. For comparison, the EMD of the quiescent 
state is also overplotted as a red, dashed line at the 
same zero level than each one of the rest EMDs. 

continuously new appearing very hot plasma 
together with heated material that evapo- 
rates from lower layers and fills the flaring 
loops. 

• During time-segments 8 and 9 the emission 
at temperatures > 28 MK turns back to the 
level found at the quiescent state, which is 
an indication of the strong heating not be- 
ing present any more. However, the height 
of the maximum in the EMD curve is even 
larger than in the previous time-segments, 
appearing also excess emission at even lower 
temperatures (down to ~ 3MK). 



4. Flare modeling 

Results from § [2] and § [3] support the idea 
of the observed rise (t > 5ks) being the con- 
sequence of two separate flare events: the first 



more impulsive one (event A) starting in time- 
segment 1, followed by a more gradual energy re- 
lease (event B) having a beginning that is merged 
with the decay phase of event A and becoming 
dominant at the end of time-segment 3. But, 
what kind of magnetic structures are pro ducing 
each one of these events? According to iReald 
([20071 ). multiple loop structures can be involved 



in a flare, but this frequently occurs only in its 
late phases. The initial phases of an X-ray flare 
are usually quite localized and one can reasonably 
assu me the presence of a single dominant loop 
(e.g. lAschwanden fe Alexander 2001 : Reale et all 



2004 ) . This assumption is realistic enough in most 
of the observed rise phases, where the impulsive 
heating typically involves a dominant loop (while 
later residual heating may be released in other 
similar adjacent loops). The clear evidence of a 
delay between the temperature peak and the den- 
sity peak is consistent with a single loop model, 
at least until the moment in which this maximum 



both in solar flares (< 
and in stellar flares (e.g 




Svlwester et al. 




993). 


. van den Oord et al.1 


1988; 


van den Oorc 


1 & Mewe 


11989: iFavata et alJ 


2000; 


Maggio et al. 


2000t Stelzer et al. 20021). even in 


very long ones (Favata et al. 20051: Getman et al. 



2008) . The presence of this delay is indeed a sig- 
nature of a relatively short heat pulse (the larger 
the delay, the shorter the heat pulse is) and of the 
coherent hydrodynamic evolution of plasma con- 
fined in a loop. Arcades and two-ribbon flares 
are instead characterized by stron g and/or last- 
ing heating ( Kopp fc Polettol 19841) and /or irreg- 
ular light curves ( jAschwanden fc Alexanderl boOl 



Reale et all 12004) . Thus, event A is likely char- 
acterized (or dominated) by a single flaring loop 
since its ignition occurs only at time-segment 1 
(as suggested by the plasma appearing at very 
high temperatures in its corresponding EMD, see 
§ 13.2. ip and shows a delay between the temper- 
ature peak and the density peak reached by the 
flaring plasma in it (see Table [2]). Event B is 
probably triggered by event A. The evolution of 
the whole EMD during event B evidences contin- 
uous, but decreasing in efficiency, heating of ma- 
terial and is consistent with event B being a two- 
ribbon flare. In § l4.2l we demonstrate that the light 
curve of event B is indeed well-reproduced by a 
two-ribbon flare model. Continual reconnections 
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in an arcade may also be responsible for event B. 
However, in this case one would expect a more 
irregular pattern in the light curve than that ob- 
served in Fig. [21 and the EMD evolution would 
not necessarily point to a kind of heating that is 
decreasing in efficiency with time, as we have ob- 
served. 

4.1. Event A: the single loop flare 

We used the results from fitting the spectra 
with the lT-model (see § [3|) , given in Table [21 to 
determine the length of the loop involved in event 
A. For the temperature in time-segment 1, we used 
kT = 4 keV (the lower limit from the fit) since its 
upper limit is undetermined (see Table [2]). 

For this event, the complete rise phase is ob- 
served and instead we miss satisfactory informa- 
tion from the decay phase. A complete theoretical 
analysis of the rise ph ase of a flare occurring in a 
single loop is given bv lReald (<2007l) . The author 
showed that the loop half-length (L) can be deter- 
mined from the maximum temperature reached by 
the flaring plasma (T ) and the temperature and 
time in which the density peaks (Tm and £m, re- 
spectively) by: 



L w 3 • 10 5/ V 



T, 



5/2 



T 2 

1 M 



(1) 



where all the parameters are given in c.g.s. units. 
The maximum temperature (To) is related to the 
temperature measured for the flaring component 
by spectral fitting (T b s )- For the EPIC instru- 
ment: 



T = 0.13T" 6 



(2) 



For event A, the maximum temperature and 
maximum emission measure were reached in the 
time-segments 1 and 2, respectively. For a single 
loop description, the (square root of the) emis- 
sion measure becomes a good proxy of the den- 
sity, because the loop volume presumably does 
not change much during the event. In this case, 
from Eq. [T] we obtain a loop half-length L = 
1.8 ± 0.3 x 10 11 cm, i.e. ~ 4 ± 1 R* (assuming 
a stellar radius i?* = 0.64 — from the pre-main 



sequence models of Si ess et al. 2000, for an M2.5 
star with 8 Myr, such as TWA 11B is supposed 
to be). Although this value is relatively large, 



Table 4: Parameters of t he flar i ng loo p involved 
in event A, derived as in iReald ([2OO70 . We have 
assumed a semi-circular geometry. 



Parameter 


Units 


Value 


n M 


10 10 cm" 3 


2.0 ±0.7 


^avg 


10 10 cm" 3 


4.1 ±0.7 


V 


10 31 cm 3 


2.3 ±0.8 


A 


10 20 cm 2 


0.6 ±0.2 


r 


10 10 cm 


0.4 ±0.1 


L 


10 11 cm 


1.8 ±0.3 



similar loop sizes have already bee n derived for 



youn g stars in star-forming r egions (jFavata et al 
2005 ; Franciosini et al.l 120071 ) . As far as we are 



concerned, this would be the first time that such a 
relatively long loop is detected in a star older than 
3 - 4 Myr. 

This loop wo uld fill a volume V ~ 2.3 x 10 31 cm 3 
fsee lReale|[2007 . for further deta Is on the relations 
used for determining this quantity). Assuming a 
semi-circular geometry, its aspect (r/L, with r be- 
ing the loop cross-section) would be of the order 
of 2 %, which is quite lower than that observed in 
th e Sun (~ 10 %), but compatible with the results 



Favata et al 



()2005l ). In Table [4] we summarize 
the main parameters of the loop involved in event 
A, i.e. those already mentioned and maximum 
density at the loop apex (tt-m), average density in 
the loop when the maximum density is reached at 
the apex (n aV g) and cross-section area of the loop 
(A). Such para meters were d etermined with the 
relations given in iRealel ( 2007 ). 



Co mparin g our results with those obtained by 



Reald (|2007Q for events observed in Algol, AB Dor, 
and Prox Cen, we conclude that: 1) our densities 
seem to be similar to those determined in the Al- 
gol and Prox Cen flares; and 2) the loop volume 
and its subtended area are of the same order of 
magnitude as those found for Prox Cen, another 
M star. Note that the largeness of some of the 
errors shown in Table [4] is a consequence of the 
propagation of errors in the equations. 

For the sake of completeness, we made the ex- 
ercise to repeat the same analysis to the entire 
flare event (A+B) as if it all occurred in a single 



9 



loop. For a flare temperature peaking in time- 
segment 1 (to ~ 6.5 ks), and the density peaking 
in time-segment 8 (<m « 37 ks) , using Eq. [IJ we 
obtained L = 1.3 X 10 13 cm, i.e. L ~ 260i?* . Ob- 
viously, such a long loop would be easily destroyed 
by the stellar rotation. Clearly, this result makes 
no sense. 

4.2. Event B: the two-ribbon flare 

For the study of event B, we used the two- 
ribbon flare model bv lKopp fc Poletto (1 984f) ex - 



tended to the stellar case ( Poletto et al. 1988) 



This model supposes that a disruptive event opens 
a loop arcade, being the open field lines then 
driven toward a radial neutral sheet (above the 
magnetic neutral line) where they reconnect at 
progressively higher altitudes. Thus, the contin- 
uous heating provided by these reconnections is 
capable of reproducing the temporal profile of the 
energy rate released during both the rise and decay 
phases of a two-ribbon flare. By analogy to that 
observed on the Sun, the model considers that the 
arcade of loops is extended along the East- West 
direction (i.e., axial symmetry around the polar 
axis is assumed). It also assumes that the mag- 
netic field is potential between the stellar surface 
and the location of the neutral line and extends 
radially outwards from there. The magnetic field 
in the meridional planes of the arcade can there- 
fore be expressed in terms of a single lobe of a 
Legendre polynomial of degree n. Note that: 

i Each lobe is latitudinally bounded by radial 
magnetic fields. 

ii The arcade corresponds to one lobe axisym- 
metrically continued over some longitude in 
the East- West direction. 

iii Through an appropriate choice of n, one can 
find a lobe placed in the range of latitudes 
covered by the active region. However, as 
spatial information is available only for the 
Sun, stellar flaring regions are generally as- 
sumed to be centered on the equator for odd 
n and to end at the equator for even n. 

iv As stellar observations cannot provide any 
information on the time-dependent rise of 
the neutral point, it is assumed to mimic 
the solar case. Thus, it follows an exponen- 
tial law of the form given by Eq. [3J where 



y is the height of the neutral point (in units 
of R+, measured from the star's center), t 
is the time (measured from the beginning 
of the two-ribbon flare, ij n i), to is a time- 
constant, and H m is the maximum height 
reached by the reconnection point during its 
upward movement (measured from the star's 
surface). H m is typically chosen to be equal 
to the latitudinal extent of the arcade, which 
is in turn linked with n (see Eqs. S]and[S]). 



I/ = 1 + ^(1- e £) 



71+1/2 



(3) 

for n > 2 (4) 
for n = 2 (5) 



Under all these assumptions, the rate of mag- 
netic energy released by the reconnecting arcade 
per radian of longitude (dE/dt) can be expressed 
as: 



dE 1 

2n{n + l){2n+lfRlBl 

dt 8ir 

1] 



/i, 2 (n) y 2n [y 



2nr,,2n+l 



P„ 2 (0i, 2 ) [n +(»+!)« 



2n+113 



1 1 - 



where Ii^{n) — J P^(0)d(cos 9) evaluated be- 
tween the latitudinal borders of the lobe, P n (6) 
is the Legendre polynomial of degree n, 6 is the 
co-latitude, is the co-latitude of either bound- 
ary of the lobe, and B m is the maximum surface 
magnetic field in the active region. The factor 
B\ x merely defines the normalization of the en- 
ergy release light curve, while to and n determine 
its shape. As Eq. |6] is given per radian of lon- 
gitude, a length (I) must be assumed for the ar- 
cade in order to calculate its total energy-release 
rate. Solar two-ribbon flares occur in loop arcades 
whose length is typically about 1.5 times their 
width (I r; 1.5H m ). In our study, we adopted this 
ratio for any given n. 

Since, in the view of many authors, the initi- 
ation process for the flare itself might be the re- 
sult of a more rapid (nearly explosive) reconnec- 
tion than the reconne ction pr ocess about which 
iKopp fc Poletto! (|l984h and IPoletto et all (|l988h 
speak, they stressed that their model is applicable 
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only after the initial flare trigger mechanism is ter- 
minated. In fact, the model seems not to be able to 
describe the time of impulsive heating and steeply 
increasing temperatures. However, the impulsive 
phase typically ceases before reaching the 50 % 
flare peak level in the soft X-ray bandpass com- 
monly used for stellar observations. At this time 
in our observations, the temperature is gradually 
decreasing. Hence, the model is applicable from 
time-segment 4 of the light curve to the end of our 
observations (time-segment 3 is rejected also for 
avoiding possible contamination from the decay of 
event A). 

In order to apply the described model, the pho- 
tospheric magnetic held in the flaring region (from 
which B m can be determinated), and its latitudi- 
nal location and size (which dictate n) should be 
known. In the solar case, observations provide all 
these data, whereas in the stellar case, at best, 
they can only be inferred indirectly. Since we are 
unaware of the location, size, and magnetic held 
strength of the active region that we are study- 
ing, we treated n and B m as free parameters to be 
determined from the best fit of the model to the 
observations. 

We created a grid of values for the free pa- 
rameters n, to, and t ml to fit event B with the 
two-ribbon flare model (the lower limit of t- ln i was 
fixed at the beginning of the whole enhancement 
because of obvious physical reasons) . For each set 
of these parameters, we determined the B m that 
best fit the data by minimizing the x 2 value, which 
is defined as: 



X 



N 

E 

i=l 



L 



mod, i 



-L. 



obs.i 



AL 



bs.i 



(7) 



where L mo d is the expected luminosity from the 
model (i mo d = / -q-l-dE/dt), L \, s is the observed 
luminosity, Ai b s is the error in the observed 
luminosity, and N the number of time-intervals 
with which we fitted the two-ribbon model (from 
time- segment 4 to 9). Fo llowing the solar anal- 



ogy (ICanfield et al.1 Il980l see § EJ), we assumed 



the measured radiative losses in the X-ray band 
to be w 10 % of the thermal energy generated 
as conseque nce of magnetic reco nnections, that is 
/ w 0.1 (see iPoletto et al.lll988l for details). Ac- 
tually, only a fraction (q < 1) of the liberated 
magnetic energy is indeed used to heat the plasma 



(thermal energy that is subsequently lost via ra- 
diation and conduction), while the rest is trans- 
formed into mechanical energy, into fast particles 
ejected from the corona, etc. At this point, we 
want to notice that we did not find a unique so- 
lution corresponding to a single set of parameters, 
but a number of solutions producing a good fit 

(x 2 ~i). 

In TableJSJ we show the two-ribbon flare param- 
eters for some of the good fits we have obtained. 
We plot these results together with the observa- 
tions in Fig. [7] Other good fits were found also for 
larger values of n (i.e., active regions with smaller 
width and, therefore, shorter loops), but note that 
a smaller active region needs higher surface mag- 
netic fields for reproducing a given energy-release 
rate. Thus, loop systems that reach larger alti- 
tudes - i.e., small values of n - may be more realis- 
tic in our case, although other loop configurations 
cannot be excluded (see § [SJ . 

Table 5: Two- ribbon flare parameters resulting 
from fitting the model to the observed energy re- 
lease rate from event B. 



Polynomial degree 


2 


3 


5 


10 


Region width [deg] 


90° 


53° 


33° 


17° 


H m [R*] 


1.57 


0.90 


0.57 


0.30 


Ll [R*] 


2.46 


1.41 


0.89 


0.47 


tin) [ks] 


7.0 


6.5 


6.0 


5.5 


to [ks] 


247 


192 


183 


174 


[G] 


440 


730 


1050 


1830 


v riso [km s -1 ] 


2.8 


2.1 


1.4 


0.8 


N e [xlO 10 cm" 3 ] 


1.1 


2.6 


5.1 


13.3 


x 2 


1.02 


1.02 


1.04 


1.04 



f L is the loop somi-lcngth determined from the maximum 
height H m assuming a semi-circular geometry. 

The values given in Table [3] for the maximum 
surface magnetic field were determined assum- 
ing that all the magnetic energy is used to heat 
the plasma that fills t he flaring loops (q = 1). 
Kopp fc Polettol (|l984h used q = 0.003 for a so- 
lar flare. On the oth er hand, dur i ng th e anal- 
ysis of a stellar flare, Giidel et al. ( 19991 ) found 
acceptable solutions with q < 0.05 for n = 2 and 
q w 0.01 — 0.02 for n — 3 — 4. For event B we found 
that small values of q (~ 0.01 — 0.02) imply very 
strong photospheric magnetic fields in the flaring 
region (B m ks 3 — 18 kG). Using q = 0.05, we ob- 
tained B m 2 kG for n = 2 and B m ~ 3 - 8 kG 
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40 60 
Time (ks) 



Fig. 7. — Best fit to the observations using the 
two-ribbon model for n = 2 (see Table [5]) . Note 
that each set of parameters in Table [5] produces 
similar results in the fitted region while having 
different slopes in the previous part and in the de- 
cay. The light curve segments used to fit the two- 
ribbon model (time-segments 4-9, see also Fig. [2]) 
is shown with a solid line. The dashed line is the 
model. We have also plotted the average values 
of the model in each time-interval during event B 
(dotted line) for a clearer comparison between the 
observed light curve and the model. 

for n = 3 - 10. 

5. The fluorescent Fe 6.4 keV line 

During the inspection of the X-ray spectrum in 
the different time-intervals, we observed a feature 
(excess emission not reproduced by the plasma 
model) close to the Fe K a line at redder wave- 
lengths. We identified this feature as the Fe fluo- 
rescent line at 6.4 keV. It was noticeable only dur- 
ing the time-segments 4 to 9, while it was neither 
in time-segments 1 to 3 nor in the quiescent state 
(see Fig. [5]). Unfortunately, we could not moni- 
tor the time-variation of this emission line since 
the individual spectra of each time-interval have 
no counts enough for accurately fitting the line. 



Table 6: Best fitting values (3T plasma model + 
Gaussian) for the spectrum integrated in the time- 
segments 4-9. Only the temperature and the 
emission measure of the third thermal component 
and the Gaussian are listed because the remain- 
ing parameters were fixed to the values estimated 
for the quiescent state. The given x^ed an d d.o.f. 
values refer to the total fit. 



kT 3 (keV) 

EMi (10 52 cm" 3 ) 

Gaussian central energy (keV) 



Gaussian flux (10 
Gaussian a (eV) 
X?ed [d.o.f.] 



ph cm 



2 02+ 013 



-0.3 

o.9±?:S 

9.7tr 7 



= 0.96 [499] 



Thus, to quantify the excess emission at 6.4 keV in 
our observations, we performed spectral fitting us- 
ing an interstellar absorbed 3T-plasma model and 
an additional Gaussian line component at 6.4 keV. 
We fitted, on the one hand, the spectrum inte- 
grated from the quiescent to time-segment 3 (in- 
clusive) and, on the other hand, the spectrum inte- 
grated from time-segment 4 to 9 (inclusive). The 
gaussian flux obtained in this way for the former 
spectrum is zero. The best fit results are shown 
in Table [5] Note that we left free only the tem- 
perature and the emission measure of the third 
(hottest) thermal component and the gaussian, 
while the rest of the parameters were fixed to the 
values obtained for the quiescent state (Table Q}. 

In the past, fluorescent Fe emission were com- 
monly o bserved in classical T Tauri stars and pro- 
tostars (|Tsuiimoto et al.ll2005t iFavata et ah 2005t 
Giardino et al.ll2007t ISciortinoll2008l ) where it has 
been attributed to the incidence of X-ray emission 
onto the proto-planetary disk or into the circum- 
stellar gas surrounding very young objects. How- 
ever, the Fe 6.4 keV line has been observed also 
in the giant star HR 9024 dTesta et al.ll2008ft and 
the RS CVn system II Peg (|Osten et alJl2007h . In 
the case of HR 9024, the authors attributed its 
presence to the incidence of hard X-rays onto the 
photosphere, in concordance with what is observed 
in the Sun. In II Peg, the excitation mechanism 
was ascribed to electron impact ionization of pho- 
tospheric Fe. 

In our observation, we have not statistics 
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Fig. 8. — Left: PN spectrum of TWA 11B integrated in the time-segments 4-9, with the clear excess 
emission in 6.4 keV. The (red) continuous line is the fitted 3T-model. In the small window, the (blue) 
dashed-line is the fitted 3T- model + Gaussian component. The Gaussian component is also plotted as a 
dotted-dashed line. Right: Same as figure at the left but in the time-segments 1-3 plus quiescent. Here, 
the excess emission in 6.4 keV is clearly not present. 



enough to perform a robust analysis of the flu- 
orescent line. Nevertheless, some constraints can 
be given. On the one hand, the line is not observed 
during the first event. This may be consistent with 
the scenario of the long loop for this event. The 
efficiency of the fluorescence decreases with the 
dist ance of the X-ray source to the photosphere 
(e.g. Drake et"al~l 2008 ). Therefore, if the fluoresce 
line were produced by photoionization, it should 
not be observed during flares occurring in long 
loops. On the other hand, the presence of the 
line during the second event could be indicative of 
a not very high loop system, if it were produced 
by photoionization. In contrast, if the line were 
produced by collisional ionization, there would 
be none constraint to the loop height. A result 
in favor with the non-photoio nization nature of 
the flu orescence line is shown in lCzesla fc Schmitt 
(J2007I). In their study, the authors modeled the 



illuminating input spectrum and obtained line 
fluxes below the observations. With our data, 
the collisional production of the fluorescence line 
cannot be discarded. 

6. Final remarks and conclusions 

In this paper, we analyzed the long rise phase of 
a flare observed in an XMM-Newton archive data 
of the - 8 Myr old star TWA 11B. The analy- 
sis of the light curve, of the hardness-ratio curve, 
and of the time-resolved spectra consistently indi- 
cates that probably the flare first involved mainly 



a single loop and then propagated to a loop ar- 
cade becoming a proper two-ribbon flare. We split 
our analysis into three parts: the quiescent state, 
the single- loop flare (event A) , and the two-ribbon 
system (event B). E vent A was st udied with the 
analysis described in iReald (l2007t). For event B 
we use d the stellar version of the IKopp fc Polettol 
(1984)'s solar two-ribbon flare model that is able 
to provide some limited information about the late 
flaring structures. 

For the single loop, we obtained a semi-length 
of approximately 1.8±0.3 x 10 11 cm (~ 4±1 R*), 
with a volume V = 2.3 x 10 31 cm 3 and a cross- 
section r/L ~ 2%. These value s are c omparable 
to those found by IFavata et all (|2005h in Orion 
members. This fact makes us suggest that large 
and thin loops are common in young active stars. 
For the two-ribbon system, different results consis- 
tent with the observed data were found. Good fits 
(x 2 ~ 1) were obtained for both small and large 
values of n (i.e. for high and short loop arcades). 
Bearing the semi-length of the first loop in mind, 
the more realistic scenario for the two-ribbon sys- 
tem is that with long loops. In any case, the es- 
timated values of the maximum surface magnetic 
field in the flaring region result to be quite strong, 
reaching 2-8 kG when a large fraction of the lib- 
erated magnetic energy (q — 0.05) is assumed to 
heat the plasma, and up to 18 kG when q ~ 0.01. 

During the inspection of the X-ray spectrum, 
we observed the Fe fluorescent line at 6.4 keV 
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during the time in which the two-ribbon system 
evolved (time-segments 4 — 9). The absence of the 
line during the first event may be consistent with 
the long loop scenario for event A. In contrast, the 
detection of the fluorescence line during the sec- 
ond event could be indicative of a not very high 
loop system involved in the two ribbon flare, if 
the line were produced by photoionization. Oth- 
erwise, there would be none constraint to the loop 
system height. 
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